Backbending in 50 Cr 
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The collective yrast band and the high spin states of the nucleus 50 Cr are studied using the 
spherical shell model and the HFB method. The two descriptions lead to nearly the same values for 
the relevant observables. A first backbending is predicted at I = 10h corresponding to a collective 
to non-collective transition. At I = 16h a second backbending occurs, associated to a configuration 
change that can also be interpreted as an spherical to triaxial transition. 



PACS number(s): 21.10.Re, 21.10.Ky 

In a recent paper [Q we have shown that large scale 
Shell Model (SM) calculations with the realistic interac- 
tion KB3 predict the same intrinsic state than mean field 
calculations (Cranked Hartree-Fock-Bogoliubov) with 
the Gogny force for the ground state rotational band 
of 48 Cr. Thus, we have two complementary views of 
the problem. For the SM wave functions have the 
proper quantum numbers (angular momentum and par- 
ticle number) and include correlations in the wave func- 
tions needed for a detailed account of the observables, 
while the mean field results provide us with a simpler 
understanding of the intrinsic state on top of which the 
rotational ground state band is built. In this paper we 
proceed in the same way for the nucleus 50 Cr. Exper- 
imentally 0, it is known that up to J = 10 the Yrast 
states of 50 Cr are members of a rotational band. The 
gamma ray energy of the last known state J = 12+ to 
the J = 10 + indicate that 50 Cr has started to backbend 
at J = 10 this situation being reminiscent of that encoun- 
tered in 48 Cr. The two additional neutrons of the nucleus 
50 Cr will contribute to the appearance of new structures 
when the angular momentum of the system is increased. 
As in the case of the Yrast band of 48 Cr, where the back- 
bending was due to a deformed to spherical transition, 
we find in 50 Cr that the incipient backbending observed 
at J = 10 is also due to a shape transition from the pro- 
late J = 10 state to the nearly spherical (actually slightly 
oblate) J = 12 state. Contrary to what happens in 48 Cr 
the Yrast band of 50 Cr undergoes a second transition at 
J = 16 from the slightly oblate regime to a more de- 
formed, triaxial regime that causes a second backbend. 
The behaviour of the Yrast band of 50 Cr can also be un- 
derstood in terms of the changes of the occupancies of 
the I/7/2 orbitals of protons and neutrons. 

Computational procedures. We shall briefly de- 
scribe the computational procedures used in this pa- 
per. A more detailed account can be found in [Q . In 
the Spherical Shell Model (SM) 50 Cr is described in a 
Ohiv space, i.e. ten particles are allowed to occupy all 
the states available in the pf shell (~ 10 7 ). The effec- 



tive interaction is the same used in the Cr calculation 
i.e., a minimally modified version of the Kuo-Brown's 
G-matrix ||] denoted KB3 in [Q. The single particle en- 
ergies are taken from the 41 Ca experimental spectrum. 
The effect of core polarization on the quadrupole proper- 
ties is taken into account by the use of effective charges 
q-x = 1.5, q v = 0.5. These reasonable values agree with 
the estimation of core effects in the quadrupole properties 
carried out for 48 Cr Q. The secular problem is solved 
using the code ANTOINE g, a very fast and efficient 
implementation of the Lanczos method. 

In the intrinsic frame calculations we have used the Self 
Consistent Cranking Hartree-Fock-Bogoliubov method 
(HFB) with the density dependent Gogny force Q. The 
mean field intrinsic states \<j>J) have been expanded in 
a triaxial harmonic oscillator basis \n x n y n z ) with differ- 
ent oscillator lengths. Ten oscillator shells are included 
in this calculation in order to ensure the convergence of 
the mean field results. In our calculation we use the DS1 
parameters set of the Gogny force J7J. Without further 
changes, this force has proven capable of describing suc- 
cessfully many phenomena, and in particular high spin 
behaviour g. 

The Yrast band. In fig. |l| the SM, HFB and experi- 
mental gamma ray energies E^(J) — E( J) — E{ J — 2) are 
plotted as a function of the angular momentum J for the 
Yrast band. The SM results are very close to the exper- 
imental data. The CHFB results follow the trend of the 
SM ones but are shifted downwards in energy indicating a 
much bigger static moment of inertia than in the SM and 
the experiment. In it was argued that this behaviour 
is the result of a deficient treatment of pairing correla- 
tions by the HFB method in the weak pairing regime. 
However, this deficiency does not affect to the nature of 
the intrinsic state in a substantial way. The experimental 
results from the Chalk- River Mc Master collaboration [El , 
show preliminary evidence for Yrast states (14+ and 16^) 
at excitation energies 9.9 MeV and 13.9 MeV. These re- 
sults would confirm the existence of a first backbending 
in 50 Cr. Notice that the corresponding E 7 's, 2.3 MeV 
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and 4.0 MeV fit reasonably well into our predictions. 

In both theoretical calculations a second backbending 
is predicted to take place at J = 16h. Its origin will be 
discussed below in terms of the changes observed in the 
quadrupole moment and of the evolution of the occupan- 
cies of the relevant spherical orbits. 

The quadrupole properties obtained in both ap- 
proaches are summarized in Fig. ^. In the upper panel 
the (3 and 7 deformation parameters obtained in the HFB 
calculation are plotted versus the angular momentum. 
We can interpret these results as defining three regions: 

- At low angular momentum 50 Cr is an axially sym- 
metric prolate nucleus (/? ~ 0.23 and 7^0) 

- At J — 10h the (3 value drops to 0.15 but the sys- 
tem is still axially symmetric while at J — 12h the 
system undergoes a shape transition to a weakly 
deformed oblate state with (3 = 0.08 and 7 = —67 
degrees. This sudden change is responsible for the 
observed backbending at J = 10ft. 

- At J = 16h a new transition occurs that makes the 
system triaxial increasing at the same time the (3 
value. These changes of the quadrupole moment 
are responsible for the backbending seen at J = 

In the middle panel the spectroscopic quadrupole mo- 
ment of the SM calculation is plotted as a function of 
J. Q s changes sign, from negative to positive at J = 10. 
This point has been recently discussed by Zamick and 
coworkers |liQ[ ]. They made truncated shell model calcu- 
lations using also the interaction KB3. Their results for 
the quadrupole moments at truncation level i=3 are qual- 
itatively equivalent to our full space results. Neverthe- 
less, they interpret these positive quadrupole moments 
as pertaining to a "high K prolate" band Jll]] while we 
rather rely in the HFB results to conclude that these 
states are of oblate, non-collective character. This tran- 
sition is correlated with the backbending in the E~ l plot. 
The spectroscopic quadrupole moment and the intrin- 
sic deformations (3 and 7 behave in a very similar way. 
However, it is difficult to establish a close connection be- 
tween them in the high spin region where the K labeling 
is probably meaningless. In ref. |Io(| it was pointed out 
that the very large and negative values of Q s for the 
states of highest spins correspond to those for the maxi- 
mally aligned configurations. Our results agree with this 
interpretation. Triaxiality comes in because two or more 
of these configurations are simultaneously present in the 
wave function of these yrast states. 

In the lower panel of Fig. || the B(E2) transition prob- 
abilities from the SM and HFB calculations are plotted 
as a function of J and compared with the available ex- 
perimental data. The HFB results are obtained from the 
intrinsic values of the quadrupole operators Q20 and Q22 



using an improved rotational formula JL2|. Both results 
are very similar in the whole range of angular momentum 
considered. The comparison with the experiment data is 
good for the B(E2, 2 + — > + ) but both theoretical results 
overestimate the experimental ones at J = 4 and 6h. 

In order to understand the behaviour of the previously 
considered observables as a function of J it is convenient 
to analyze the structure of the intrinsic states in terms of 
the "fractional shell occupancies" v(n, I) and the "shell 
contribution to (Jx)" j x (n : l) defined in ref. jlj]. In Fig- 
ure [| we have plotted the fractional occupancies of the 
spherical orbits in the HFB solution (upper panel) and 
in the SM one (middle panel). In addition the shell con- 
tributions to (J x ) (which only make sense in the HFB 
calculation) are plotted for the relevant spherical orbits 
in the lower panel. There are striking similarities between 
the SM and HFB occupancies, with the HFB results be- 
ing slightly smoother than the SM ones. Let us examine 
these numbers for the three regions sketched before. 

Up to J = 8h only the I/7/2 orbit and its Aj = 2 de- 
formation partner 2p 3 / 2 are appreciably occupied. This 
is in agreement with the fact that up to that value of the 
angular momentum the nucleus is an axially symmetric 
prolate rotor {(3 ~ 0.23) as surmised in ref. [jl3). The 
situation is the same than in 48 Cr. For J = 8fr the oc- 
cupation of the 2p 3 / 2 orbit starts decreasing in favor of 
the I/7/2, as the later is more capable to create angular 
momentum. As a consequence the (3 deformation gets 
reduced (see Fig. |^) and the first backbending sets in. 

At J = 12h the 2p 3 / 2 orbit is emptied both for pro- 
tons and neutrons eliminating the quadrupole coherence. 
In the lower panel we can see that the contribution of 
the neutrons to the total angular momentum, J = 12h, 
is J — 6h and comes fully from the I/7/2 orbit. But, 
six units of angular momentum is the maximum attain- 
able by six identical particles in the I/7/2 orbit, therefore 
neutrons are fully aligned. In order to build up fourteen 
units of angular momentum the system chooses to align 
the four protons in the I/7/2 orbit to obtain J = 8h (the 
other alternative is to move neutrons to higher orbits). 
J = 14h is the maximum angular momentum that can 
obtained with four protons and six neutrons in the 1/7/2 
orbit. Therefore, to reach J = 16h, the system is forced 
to pump out particles of the low K states of the I/7/2 
orbit to high K states of the I/5/2 as seen in the neutron 
occupancies plot. This leads to the second backbending 
because the excitation energy of the 16 + state is anoma- 
lously large. When the I/5/2 neutron shell is exhausted 
(at J = 187i) the 2p 3 / 2 neutron shell enters into the game 
again. Placing particles in the 2p 3 / 2 increases again the 
quadrupole coherence producing the increase of deforma- 
tion observed in Fig. || At J = 22?i the band terminates 
as the maximum angular momentum that can be built 
with four protons and six neutrons in the fp shell is 23?i. 

In figure ^ we present our predictions for the gyromag- 
netic factors compared with the experimental results of 
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Pakou et al. |14j . Once again the SM and HFB results are 
nearly identical for the whole band but they do not agree 
with the data except for the state J = 2, The situation 
is somehow puzzling. We know that 50 Cr is not as good 
a rotor as 48 Cr; its B(E2)'s are smaller and the collec- 
tivity, as a function of the rotational frequency vanishes 
earlier. Our descriptions of 50 Cr tend to overshoot the 
experimental B(E2) values, indicating that possibly we 
are obtaining a too large quadrupole collectivity. How- 
ever, for the gyromagnetic factors, the contrary seems 
to happen. We need to increase the quadrupole correla- 
tions in order to approach the rotational limit indicated 
by the experiment. It is difficult to move simultaneously 
in both directions. We have tried another version of the 
KB interaction called KB' in jlB) that has a smaller gap 
between the 2p 3 / 2 and the I/7/2 orbits, thus favoring de- 
formation. As we feared, while the g values get reduced 
(c?(2+)=0.53, 3 (4+)=0.60, £r(6+)=0.62, 5 (8+)=0.66) and 
nearly agree with the experimental results of ref . ||14| , the 
B(E2)'s increase a lot. For instance the 4+ — > 2 + transi- 
tion has B(E2)=415 e 2 fm 4 , compared to the experimen- 
tal value 159(21) e 2 fm 4 or to the KB3 value 264 e 2 fm 4 . 
Therefore we think that it will be extremely difficult to 
reproduce simultaneously the experimental g factors and 
the experimental B(E2)'s. 

Further discussion of the high spin region. We 

have already shown that the double backbending in 50 Cr 
is related to the existence of three different zones in the 
yrast line, and takes place in the two interphases. In the 
first of those regions, up to J = 10H the decay proceeds 
by the usual sequence of A J = 2 enhanced E2 transitions. 
Here we shall study the decay patterns in the other two 
zones. 

In figure ^| we have plotted the yrast band as given 
by the SM calculation, from the first backbending to the 
maximum J attainable in the (I/7/2) 10 configuration. A 
new feature is evident: beyond J — 10h the yrast se- 
quence shows up as A J = 1. Therefore, the decay may 
become more complex, in particular Ml transitions are 
now allowed. And not only allowed but dominant due 
to the hampering of the quadrupole collectivity in this 
zone. In the figure we can follow the decay subsequent 
to the feeding of the J = 14?i state. Indeed the emitted 
gammas are mostly Ml, only the 14+ — > 12 + is domi- 
nantly E2. The figure also presents our predictions for 
the lifetimes and branching ratios using the measured 
E 7 's. Experimental results from the Mc Master Chalk 
River collaboration and from Legnaro fl(i}| are expected 
to come out soon. 

It is quite exciting to find out that the yrast decay at 
the first backbending proceeds via magnetic transitions. 

Let's now move higher in energy and angular momen- 
tum. If we assume that the yrast band is fed at the yrast 
18+ state, the dominant decay sequence will be: 
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If we choose the yrast 17"* 
decay sequence will be: 
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Hence, in the third zone we recover the AJ = 2, E2 
decay pattern, in spite of the fact that the yrast sequence 
is now 14+ - 15+ - 16+ - 17+ - 18+. The surprising thing 
is that the decay bypasses the yrast 14+ and 13+ states, 
a striking feature associated to the second backbending. 

In conclusion, we have shown that the spherical 
shell model and the Cranked Hartree Fock Bogolyubov 
method give similar -and accurate- descriptions of 50 Cr. 
In particular, they predict two backbendings in the E 7 
plot of the yrast band, which are interpreted as due 
to two changes of regime: first from collective pro- 
late to non-collective oblate and second from (I/7/2) 10 - 
dominated configurations to aligned states with less than 
ten particles in the I/7/2 shell. A puzzling situation is 
found in the collective, low energy part of the yrast band, 
where the calculations fail to reproduce the experimental 
gyromagnetic factors. Finally in our study of the decay 
sequence from the high spin yrast states we have found 
out that the decay mode varies from one zone to another. 
In the high spin regime it proceeds by stretched E2 tran- 
sitions, at the first backbending by Mi's and in the lower 
part by collective E2's. At the second backbending the 
decay sequence bypasses the yrast line . 
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FIG. 1. Theoretical (triangles, SM; circles HFB) and ex- 
perimental (squares) gamma ray energies versus the angular 
momentum J 
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FIG. 2. Upper panel: j3 and 7 deformations in the HFB 
calculation. Middel panel: Spectroscopic quadrupole moment 
Q 3 computed in the SM approach. The experimental value 
Q s (2 + ) is also given. Lower panel: B(E2, J — > J— 2 transition 
probabilites versus J computed in SM (triangles) and HFB 
(circles) compared to the experimental data (squares) 
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FIG. 3. Upper panel: "Fractional shell occupancies" v(n, I) 
computed in the HFB approach as a function of J. Middel 
panel: Same as before but for the SM. Lower panel: "Shell 
contribution to (J x )" j x (n,l) in HFB. 
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FIG. 4. Gyromagnetic factors: SM (triangles), HFB (cir- 
cles) and experiment (squares). 
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FIG. 5. Decay scheme in the region of the first backbend- 
ing. Energies in keV. In parenthesis the experimental values. 
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